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Abstract The numerical generation of random quan¬ 
tum states (RQS) is an important procedure for in¬ 
vestigations in quantum information science. Here we 
review some methods that may be used for performing 
that task. We start by presenting a simple procedure 
for generating random state vectors, for which the main 
tool is the random sampling of unbiased discrete prob¬ 
ability distributions (DPD). Afterwards the creation of 
random density matrices is addressed. In this context 
we first present the standard method, which consists in 
using the spectral decomposition of a quantum state for 
getting RQS from random DPDs and random unitary 
matrices. In the sequence the Bloch vector parametriza- 
tion method is described. This approach, despite being 
useful in several instances, is not in general convenient 
for RQS generation. In the last part of the article we 
regard the overparametrized method (0PM) and the 
related Ginibre and Bures techniques. The 0PM can 
be used to create random positive semidefinite matrices 
with unit trace from randomly produced general com¬ 
plex matrices in a simple way that is friendly for numer¬ 
ical implementations. We consider a physically relevant 
issue related to the possible domains that may be used 
for the real and imaginary parts of the elements of such 
general complex matrices. Subsequently a too fast con¬ 
centration of measure in the quantum state space that 
appears in this parametrization is noticed. 
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1 Introduction 

About three decades ago Paul Benioff m and Richard 
Feynman EH envisaged a computer whose basic con¬ 
stituents could be in a complex quantum superposition 
state. In the last few years we have been witnessing 
astonishing theoretical and experimental developments 
in quantum computing and quantum simulation [51161 
[7] , and also in others sub-areas of quantum information 
science Ennnniin], with experimental implementations 
already going beyond the best present classical capabil¬ 
ities [T7] . These are the first sights of what will turn out 
to be a revolution in our science and technology [Hill]. 

Nevertheless, before that can in fact become a real¬ 
ity, we still have much to understand concerning quan¬ 
tum systems with many degrees of freedom. One im¬ 
portant tool for accomplishing this task is the gener¬ 
ation and analysis of RQS [iniiTiiiiTKiiKiaiioiiiiiiia 
l23U241l25ll261l27l[2^[29l[80l[3T| , which will have an analo¬ 
gous role to that that random numbers have in classi¬ 
cal stochastic theories The parametriza¬ 

tion of quantum states [miMiEi] is the initial step to¬ 
wards generating them numerically and is one of the 
main topics of this survey, which is organized in the 
following manner. In Sec. [5] we consider the genera¬ 
tion of random pure states, for which the availability 
of unbiased random discrete probability distributions 
is indispensable and is hence also recapitulated. The 
remainder of the article is dedicated to the creation of 
general random density matrices. In Sec. [3] the stan¬ 
dard method is described together with the Hurwitz 










2 


Jonas Maziero 


parametrization for unitary matrices, which used in its 
implementation. Subsequently, in Sec. 01 the Bloch vec¬ 
tor parametrization, though impractical for RQS gen¬ 
eration, is regarded for completeness. The last part of 
the article, the Sec. [5l is dedicated to present and in¬ 
vestigate some issues regarding the overparametrized 
and related methods. In Sec. EH we discuss unwanted 
physical consequences of the ranges usually used in the 
literature for the real and imaginary parts of the ele¬ 
ments of the general complex matrices involved in this 
method and present a simple solution for the problem. 
In Sec. 15.21 we report an important possible drawback 
of the 0PM regarding its use for random sampling in 
the quantum state space: its too fast concentration of 
measure. We discuss the Ginibre and Bures methods in 
Sec. 15.31 A brief summary of the article is presented in 
Sec. El 

2 Pure States 

When there is no classical uncertainty about the state 
of a quantum system, it is represented by a vector in a 
Hilbert space %. For discrete systems, % is simply 
with the inner product between any two of its vectors 
defined as (V’|0) := where z* is 

the complex conjugate of z and d is the system dimen¬ 
sion. Here we use the standard notation of Dirac for 
vectors and, for a generic matrix A, we denote as its 
adjoint (conjugate transpose). Any state \tj;) G H can be 
written as a linear combination of the vectors of any ba¬ 
sis. One base of special interest is the computational or 
standard basis: |ci) = [10 • • • 0]^, [ 02 ) = [0 1 • • • 0[^, 
• • •, \cd) = [0 0 • • • 1[^, in terms of which 

W = EU^j\cj), (1) 

with i/jj = (cjl'tjj)- Above, X'^ denotes the transpose of 
the matrix X. 

The Born’s probabilistic interpretation of the state 
vector \tjj) requires its normalization: 

\m--=Vm) = \/TU\^=^- (2) 

Thus, as the numbers are non-negative and sum 
up to one, they form a probability distribution pj := 
Using ipj = lipjl expiiOj) one can write 

W = exp(i6»j)[cj), (3) 

with the phases 6j G [0, 27r]. 

Now we recall that if we have access to a random 
number generator yielding random numbers with uni¬ 
form distribution in [0,1], an unbiased random discrete 


probability distribution (RDPD) pITlIIIH] can be gener¬ 
ated as follows m- First we create a biased RDPD gen¬ 
erating qi in the interval [0,1] and qj in [0,1 — Pk] 
for j = 2, • • • , d. Then we use a random permutation 
of {1, • • • , d}, let us call it {fci, • • • , kd}, and define the 
unbiased RDPD as 

{Pi,--- ,Pd} := {dfei,--- ,dfcd}- (4) 

The unbiased RDPD generated in this way and d 
independent random phases 6j are then applied to gen¬ 
erate a random pure state. It is worth observing that 
there will be no privileged direction in "H only because 
the RDPD is unbiased. This pure state generation pro¬ 
cedure gives \ip) distributed with a Haar measure. An¬ 
other manner of obtaining samples with similar proper¬ 
ties is by using the rows or columns of random unitary 
matrices, which we shall discuss in the next section. 

3 Standard method 

The states of a d—level quantum system are described, 
in the most general scenario, by a density matrix p |401 
HU, which is a Hermitian positive semidefinite matrix 
(notation: p > 0) with unit trace (Tr(p) = 1). Any 
such matrix can be written in the form of a spectral 
decomposition: 

P = E,tiGlG)(Gl, (5) 

with the real eigenvalues of p being nonnegative (rj > 0 
for all j = 1, • • • , d) and summing up to one (X]j=i G = 
1). That is to say, is a probability distribution 

[SUES]. The eigenvectors of p, {\rj)}'j^n form an or¬ 
thonormal basis for the vector space i.e., {rj\rk) = 
Sjk and J2j=i Ig)(g I ~ where is the dxd identity 
matrix. 

Let us briefly look over the number of real param¬ 
eters needed for a complete description of an arbitrary 
density matrix. In order to describe the probability dis¬ 
tribution the eigenvalues of p, we need d — 1 

real numbers. Besides, as any two bases for the vector 
space are connected by an unitary matrix U (i.e., 
UC/t = Id), one can write 

Ig) = U\cj), (6) 

for j = I,-- - ,d, with {[cj)}^^]^ being the computa¬ 
tional basis, as shown in Sec. 0] Therefore the bases 
{lG)}i=i completely determined by U. Once d? — d 
real parameters are sufficient to specify completely an 
arbitrary unitary matrix U with dimensions dxd [36], 
it follows that d^ — 1 independent real parameters are 
sufficient for a thorough description of any density ma¬ 
trix. 
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From the last two paragraphs, we see that the nu¬ 
merical generation of a RQS (using the density matrix 
as written in Eq. ®) can be cast in terms of the cre¬ 
ation of a RDPD and of a random unitary matrix (RU) 
[36ll^l4^l44] . We call it the standard method because it 
would be a natural first choice giving the defining prop¬ 
erties of a density matrix. Moreover, it utilizes as few 
real parameters as possible. This is a nice characteris¬ 
tic in the view that for doing some statistics with RQS, 
one in general needs to generate many of them, what 
can be a very time-consuming task for large values of 
the system dimension d. 

From the several possibilities available [30] > in this 
article we choose the Hurwitz parametrization for gen¬ 
erating RUs [43] . In this parametrization one writes any 
dx.d unitary matrix U in terms of unitaries 

( 7 ) 

in bi-dimensional sub-spaces. The non-null elements of 
such elementary transformations are: 

^k'k^ = 1 for fc = 1 , • • • ,d and k ^ i,j] 

= cos((/)ij)exp(fV'y); U^y'’ = sin((()y)exp(fxy); 

Uff = cos((?iij)exp(-iV'ij)- ( 8 ) 

A general unitary transformation, for a d-level quan¬ 
tum system, can then be written as 

U = eMia)UiU2U^---Ud-u ( 9 ) 


with the sub-matrices being 
Ui = ' 012 , X 12 ), 

U 2 = 023, 013, X 13 ), 

c /3 = C/(3d)(<^34,i/,34,0)C/(2d) (024, 024,0) 

C/d_i = C/(‘^-i’‘^)(0d_i,d,0d_i,d,O) 

jj{d-2,d) 0d-2,d, O)--- 

t/('’"H 01 ,d, 01 ,d,Xl,d). ( 10 ) 

The random numbers appearing in the last equations 
are distributed uniformly in the following ranges of val¬ 
ues: 


0 < a < 27r; 0 < ipij < 27r; 0 < Xij < 2'?r; (H) 

(j)ij = arcsin(^V^*), 0 < < 1, f = 1, 2, • • • , d - 1. 

It is worthwhile mentioning that, although not ad¬ 
vantageous, it is possible to use the rows or columns of 
such a random unitary matrix as random state vector. 


4 Bloch vector parametrization method 


The Hermitian-traceless-orthonormal generators of the 
special unitary group SU{d), Fj {j = 1, - ■ ■ — 1), 

and Id can be used as a basis in terms of which we 
can write any dxd density matrix in the so called Bloch 
vector parametrization |36j : 

d^-i 

( 12 ) 

j=o 

where 7 = ( 71 , • • • , 7 ^ 2 - 1 ) is the so called Bloch’s vec¬ 
tor. One can use Tr(p) = 1 to see that 70 = 1/d and 
Tr(FjFk) = 2Sjk to show that the coefficients in Eq. 
m are half of the mean values of the aforementioned 
generators of SU{d), i.e., = 2“^(Fj) G K. 

For producing random quantum states using the 
Bloch vector parametrization, d^ — 1 real random num¬ 
bers 7 j must be generated. The main difficult here is 
that for d > 3 there is no known explicit determination 
of the range of values for the parameters Xj that will 
lead to a physical state. Thus, given a basis for SU{d), 
we may use the spectrum of each Fj to determine the 
range from which we shall sample the corresponding 
7 j. In the context of RQS generation, one attractive 
choice for the generators of SU{d) are the generalized 
Cell Mann matrices: 


|ci)(cfc| -b \ck){cj\ for 1 < j < fc < d, 
-i\cj){ck\ +i\ck){cj\ for 1 < / < fc < d, 

\cj){cj\ - /|q+i)(q+i| ^ ^ ^ 

— - - , - for 1 < Z < 

vw+m ~ ~ 


1 . 


(13) 

(14) 

(15) 


A simple analysis shows that for the generators in Eqs. 
m and m we have xj € [— 1 / 2 , 1 / 2 ] while for those 
in Eq. m I 3 e [-VV(2(^ + 1)), l/vW + ^j- 

Although the condition Tr(p) = 1 is promptly satis¬ 
fied, after generating the entire Bloch’s vector we must 
yet do a positivity test. This task requires much com¬ 
putational time, what makes this method impractical 
for the task under scrutiny here. 


5 Overparametrized method 


The basic motivational idea for this method comes from 
the simple observation that, for any complex matrix 
A = (Ajfc), we have: (0|AlA|0) = ||A|0)||2 = |||0 )||^ > 
0 , where | 0 ) is any vector of and ||| 0 )|| := -y/ { 4 >\ 4 >) 
is the Euclidean norm of the vector |0) G C'’*. That is 
to say, for a general complex matrix A, the matrix A'^A 
is guaranteed to be positive semidefinite (A^A > 0 ). 
Thus, if A is normalized, i.e., if we define 



(16) 
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it is possible to write a valid density operator as: 

p = A^A. (17) 

Above ||A ||2 := a/(A| A) is the Hilbert-Schmidt norm of 
A, with {A\B) := Tr(Ali3) being the Hilbert-Schmidt 
inner product between the matrices A and B [41]. 

The simple formula for p in Eq. m has found ap¬ 
plications in quantum information science |151l2dl[!?71l451 
[46|. Once the complete description of a general com¬ 
plex dxd matrix A requires real parameters, one 
notes that this parametrization, despite being simple 
and friendly for numerical implementations, uses more 
real numbers than necessary, as discussed above. Thus 
it is dubbed as the over parametrized method. The nu¬ 
merical generation of RQS via this method is further 
explained in the next sub-section. 

5.1 An issue on the domains of Re(Ajfc) and Im(Ajfe) 

Let us start our analysis of the production of RQS via 
the overparametrized method by considering the sim¬ 
plest quantum system, a two-level system also known as 
quantum bit, or qubit for short. The advantage of using 
this system as our starting point is that it can be visu¬ 
alized straightforwardly in the . For that purpose we 
simply need to write a density operator p using the 2x2 
identity matrix I 2 and the Pauli matrices (Jj (j = 1,2,3) 
as a basis (the case d = 2 in the Bloch method): 

p = 2 -A 2 + (18) 

where xj = Tr(p(Tj) is the value of the component of 
the system’s “polarization” in the direction j = 1, 2, 3 = 
X, y, z. The real numbers (a;i, X 2 , X 3 ) = {x, y, z) are used 
as the Cartesian coordinates in R.^. Enforcing the p in 
Eq. (IT^ to be a density matrix leads to the following 
restrictions [30]: — 1 < < 1 and X]j=i x'j < 1. There¬ 

fore the points (xi, X 2 , X 3 ) must lie within a ball with 
radius equal to one and centered at (0,0,0), known as 
the Bloch’s ball (BB). 

There are several functions one may be interested in 
when working in quantum information science. Some 
relevant examples are quantifiers for total correlation 
m, quantum entanglement |3H] j quantum discord [31] , 
quantum coherence [sniEii, and quantum channel ca¬ 
pacities m- All of these quantities can, in general, be 
defined using distance measures in the quantum state 
space. For our purposes in this article, the Hilbert- 
Schmidt distance (HSD) fits well. The HSD between 
two density matrices p and is defined as the Hilbert- 
Schmidt norm of their subtraction 133111]: 

(19) 


If Xj are the real eigenvalues of the Hermitian matrix 
p — then 

dnsip, C) = v/Tr((p-C)np-C)) = (20) 

For the calculations involved in this article, the so 
called Mersenne Twister method [S3] is used as the 
pseudo-random number generator (pRNG) and the LA- 
PACK subroutines [Sj are utilized for computing eigen¬ 
values. With these tools at hand, when the standard 
method described in Sec.[3]is applied for generating one- 
qubit pseudo-random quantum states, the distribution 
of such states in the Bloch’s ball and the histogram for 
the probability of the possible values of HSD are those 
shown in the upper green panel of Fig. [T] It is worth 
mentioning that the higher density of states observed 
closer to the center of this figure can be understood by 
noticing that the direction in R^ defined by U is random 
and that ri and r 2 = 1 — ri are uniformly distributed 
in the interval [0,1]. 

Let us consider the same kind of computation, but 
applying now the overparametrized method for generat¬ 
ing the pseudo-RQS. For that purpose the pRNG can be 
utilized for the sake of obtaining pseudo-random num¬ 
bers for generating the real, 

Re(A,fc) =: (21) 

and imaginary, 

Im(A,fe) =: AJfe, (22) 

parts of the matrix elements of A = (Ajk). The first 
issue we want to deal with here is with regard to the 
domains that one may use for those numbers. For in¬ 
stance, we can follow Refs. [HIISZIISS] and generate the 
matrix elements Ajk using uniformly distributed ran¬ 
dom numbers and setting 

A”„ AJ.e [0,1]. (23) 

As shown at the right hand side of the gray panel 
at the middle of Fig. [T] the probability distribution 
for the HSD obtained in this way is, to some extent, 
qualitatively similar to that obtained using the stan¬ 
dard method. This may lead to the impression that 
our choice for the domain of the matrix elements is 
fine. However, a rapid inspection of the distribution of 
states in the Bloch’s ball obtained using the 0PM with 
A”^, Aji. S [0,1] reveals a misfortune. Even though the 
polarization in the y and 2 ; directions have approxi¬ 
mately equal chance to be positive or negative, only 
positive values for the polarization in the x direction 
are generated. There is no need to say that such a re¬ 
striction over the possible values of physical observables 


dhs{pX) ■= l|p- CIb- 
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Fig. 1 (color online) One the left is presented the distri¬ 
bution in the Bloch’s ball of two thousand pseudo-random 
one-qubit states generated using standard method (upper 
green panel) and using the overparametrized method with 
the ranges for the matrices elements as utilized in Refs. [T51 
1371 . i.e., A^f., G [0,1] (gray panel in the middle) and 
with A^i^, E [—1,1] (pink panel at the bottom). One 
the right hand side is shown the probability distribution for 
the Hilbert-Schmidt distance of one million pseudo-random 
quantum states generated using the corresponding method 
(see the text for more details). 


of the system is not a desirable feature for a method 
supposed to generate random quantum states. 

We notice that a simple solution for this problem is 
generating the matrix elements Ajk with 

(24) 

With this change, for this case, the distribution of states 
in the BB becomes even more uniform than that that 
we get using the standard method, as shown in the pink 
panel at the bottom of Fig. [TJ We want to emphasize 
already at this point that increasing the range of values 
for A^f, and does not causes any significant modi- 
hcation neither of these results nor of those that shall 
be reported in the next sub-section. 


5.2 A too fast concentration of measure for the 0PM 

In the previous sub-section we showed that the applica¬ 
tion of the overparametrized method with the real and 
imaginary parts of Ajk drawn randomly and uniformly 
from the interval [—1,1] yields an uniform distribution 


of one-qubit pseudo-random density matrices. This en¬ 
couraging result leads naturally to the question of if 
such a scheme can be applied appropriately for random 
sampling in high-dimensional quantum systems. In this 
section we investigate this question and present strong 
evidences for answering it in the negative. 

It is known for some time now that in high di¬ 
mensional spaces random variables tend to concentrate 
around their mean values |56] . In the last few years, 
this phenomenon of concentration of measure, that is 
formalized in Levy’s Lemma, has gained great impor¬ 
tance and utility in quantum information science (see 
for instance Ref. m and the references therein). 

Notwithstanding, as shown in the gray panel on the 
right hand side of Fig. [51 the 0PM leads to a too fast 
concentration of measure for the Hilbert-Schmidt dis¬ 
tance in the quantum state space as the system’s di¬ 
mension d increases. We note that such a concentra- 
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Fig. 2 (color online) Probability distribution for the differ¬ 
ent possible values of the Hilbert-Schmidt distance for one 
million pairs of quantum states generated using the standard 
method (green panel on the left) or generated via the over¬ 
parametrized method with Aj^., Aj^. E [—1, 1] (gray panel on 
the right). We see that, in contrast to what happens in the 
standard method, there is a too fast concentration of mea¬ 
sure in the OPM as the system’s dimension d increases. We 
observe that although only some values of d are shown in this 
figure (1, 2, 3, and 4 qubits), the mentioned effect is smooth 
and gradual. 
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tion of measure is much more slow in our benchmark¬ 
ing method: the standard method with the Hurwitz’s 
parametrization for unitary matrices. It is worth ob¬ 
serving that, in part, the shift in the probability distri¬ 
bution for the HSD observed with the standard method 
(green panel on the left hand side of Fig. [5]) can be un¬ 
derstood as being due to the fact that as d increases 
the same number of points will be spread in a “bigger” 
space, diminishing thus the chance for closer pairs of 
configurations to be generated. 

It is important mentioning that if instead of gener¬ 
ating A as described above, we draw it from the Ginibre 
ensemble [58ll5^ . i.e., if we produce >1”^ and using 
random numbers normally distributed (and with aver¬ 
age equal to zero and variance equal to one), very sim¬ 
ilar results are obtained, as is shown in Table [TJ Thus 
the effect seems to be a characteristic trait of the over¬ 
parametrized method, being independent on how it is 
applied. 



Fig. 3 (color online) Average of the Hilbert-Schmidt dis¬ 
tance, (dha), and the associated standard deviation, Ad^a, as 
a function of the Ginibre matrix left dimension d' for some 
values of its right dimension d. A sample with one million 
pairs of dxd density matrices was created using the Ginibre 
method for each pair {d', d). We see that both (d^a) and Ad^a 
decrease with d' for a specified value of d. 



Uniform 

Normal 

Standard 

d 

{dhs) 

^dhs 

{dhs) 

^dhs 

{dhs) 

^dhs 

2 

0.697 

0.267 

0.728 

0.267 

0.524 

0.243 

4 

0.626 

0.111 

0.655 

0.113 

0.702 

0.213 

6 

0.538 

0.063 

0.558 

0.065 

0.794 

0.204 

8 

0.476 

0.042 

0.490 

0.043 

0.844 

0.195 

10 

0.431 

0.034 

0.442 

0.031 

0.874 

0.190 

12 

0.396 

0.023 

0.405 

0.024 

0.894 

0.185 

14 

0.369 

0.019 

0.376 

0.019 

0.908 

0.182 

16 

0.346 

0.015 

0.352 

0.016 

0.918 

0.179 


Table 1 Mean value ((d^s)) and standard deviation (Adha) 
of the Hilbert-Schmidt distance for one million pairs of 
d—dimensional quantum states randomly generated using the 
overparametrized method with uniformly or normally dis¬ 
tributed random numbers or generated using the standard 
method. 

We also see in Table [T] that, even though the concen¬ 
tration of measure is ubiquitous, while the width of the 
probability distribution for the HSD obtained via the 
0PM applied to four qubits is less than 6 % of that ob¬ 
tained in the one-qubit case, for the standard method 
the corresponding percentage is almost 74 %. We no¬ 
tice another bold difference between the two methods: 
as d increases, they shift (dhs), the “typical” value of 
the HSD, in opposite directions (see also Fig. [ 21 ). 

5.3 Ginibre and Bures methods 

For completeness, in this sub-section we briefly describe 
two other methods for RQS generation whose starting 
point is also the sampling of matrices from the Gini¬ 
bre ensemble. Let us begin with a generalization of the 
0PM, that will be named here as the Ginibre method. 


If the d'yid Ginibre matrix H is a square matrix as the 
ones considered in Sec. 15.21 i.e., if d' = d, the RQS 
are generated with a Hilbert-Schmidt measure. On the 
other hand, in the general case where the number of 
lines and columns of A need not to coincide, the RQS 
are said to be generated with an induced measure [S5] . 
Using the Ginibre method to generate a sample with 
one million pairs of states for each pair (d', d), we show 
in Fig. [3] the dependence with d' of the average and 
standard deviation of the Hilbert-Schmidt distance for 
some values of d. We see a strong dependence of both 
quantities with d'. This raises an additional practical 
question about this method. Which value of the Gini¬ 
bre matrix left dimension d' should be used and how to 
justify the choice? 

Now we describe the other method, which shall be 
dubbed as the Bures’ method because it leads to RQS 
with a Bures measure. This is accomplished by defining 

m 

(I, + I/)HHt(I, + Ut) 

^ Tr((Id + C/)HHt(Irf + Ut))’ ^ ^ 

with A being a dxd Ginibre matrix and U is dxd a 
random unitary matrix. It is note worthy that 3d^ — d 
real parameters are necessary to create a RQS via this 
method. For one million pairs of states generated in 
this way, we show in Fig. 0] the center and width of 
the probability distribution for the Hilbert-Schmidt dis¬ 
tance as a function of the system dimension d. A be¬ 
havior similar to that observed for the overparametrized 
and Ginibre methods, discussed respectively in the last 
sub-section and in the last paragraph, is seem here. 
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Fig. 4 (color online) Average Hilbert-Schmidt distance (d^s) 
(black points) as a function of the system dimension d. A 
sample with one million pairs of density matrices was gen¬ 
erated, for each value of d, using the Bures’ method. In the 
shadowed cyan area are shown values of the HSD standing 
between {d^s) - ^dhs and {dus) + ^dhs- 



However the rate of concentration of measure is a lit¬ 
tle less pronounced when compared with that for the 
0PM. For the Bures’ method the width of the proba¬ 
bility distribution for four qubits is approximately 8 % 
of that for one qubit. 


6 Final remarks 

In this article we presented a brief survey of some meth¬ 
ods that may be used for the numerical generation of 
random quantum states. We gave particular emphasis 
to the overparametrized method, which is frequently 
used in quantum information science. After utilizing a 
qubit system to identify and solve a physically relevant 
problem related to the domains of the matrix elements 
used so far in the literature in implementations of the 
0PM, we considered its possible application for random 
sampling in high-dimensional quantum systems. In this 
last scenario we showed that the overparametrized and 
related methods lead to a too rapid concentration of 
measure that may prevent any fair random sampling of 
quantum states, even for quantum systems with mod¬ 
erate dimension. 
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